Frustration and time reversal symmetry breaking for Fermi and Bose-Fermi systems 
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The modulation of an optical lattice potential that breaks time-reversal symmetry enables the 
realization of complex tunneling amplitudes in the corresponding tight-binding model. For a super- 
fluid Fermi gas in a triangular lattice potential with complex tunnelings the pairing function acquires 
a complex phase, so the frustrated magnetism of fermions can be realized. Bose-Fermi mixture of 
bosonic molecules and unbound fermions in the lattice shows also an interesting behavior. Due to 
boson-fermion coupling, the fermions become slaved by the bosons and the corresponding pairing 
function takes the complex phase determined by bosons. In the presence of bosons the Fermi system 
can reveal both gap and gapless superfluidity. 
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I. INTRODUCTION 

Cold atoms in optical lattices provide a unique medium 
for mimicking effects known from otlier areas of physics. 
This is primarily due to a great flexibihty and a precise 
manipulation of the cold atomic system [iHll- Atoms 
of a fermionic or bosonic character may be placed in 
an optical lattice potential whose geometry may be eas- 
ily controlled by changing directions and/or polariza- 
tions of laser beams. Interactions between atoms may 
be controlled via magnetic, optical or microwave Fesh- 
bach resonances 0, [al- The change in the depth of the 
optical lattice modifies primarily the tunneling between 
lattice sites (having also effect on the effective interac- 
tion strength) enabling e.g. the superfluid-Mott insula- 
tor quantum phase transition in the optical realization 
of the Bose-Hubbard model as proposed by Jaksch and 
ZoUer @ and subsequently demonstrated in Ref. Q- 

Another spectacular way of controlling the tunneling 
has been proposed by Eckardt, Weiss and Holthaus Q. 
Fast periodic modulations of the optical lattice allow for 
an effective, time-averaged tunneling to be totally sup- 
pressed keeping the depth of the lattice potential un- 
changed. By varying the strength of the modulation one 
can induce the superfluid-Mott insulator quantum phase 
transition Q as verified experimentally a few years later 

[loj. Importantly, not only the magnitude, but also 
the sign of the tunneling amplitudes can be altered us- 
ing this approach. This concept has been utilized in a 
recent proposition to create frustrated magnetism with 
cold bosons in a triangular lattice fll| later implemented 
in fascinating experiments of the Hamburg group [l^ . 

The effective tunneling caused by periodic lattice mod- 
ulations is adequately explained in the framework of the 
Floquet theory Q for the periodically time-dependent 
Hamiltonians. The properties of the corresponding quasi- 
energies spectra are known to depend on the global sym- 
metries of the Hamiltonian |13l - tl5fl . Employing the simi- 
lar idea we demonstrate that periodic perturbations that 
break the time-reversal invariance (TRI) can change not 



only the sign of the tunneling amplitudes, but may also 
induce their complex values. 

In the following, we concentrate on superfluid fermions 
in the Bardeen-Cooper-Schriffer (BCS) regime with bro- 
ken TRI. We show that typically a pairing function for 
s-wave interactions acquires a complex phase which may 
be controlled by the TRI breaking mechanism consid- 
ered in thepresent paper (for a discussion of p-wave or- 
bitals see [13 )■ We also point out that in the Bose-Fermi 
mixture with broken TRI the complex ground state of 
bosons affects the Fermi pairing function. The effect is 
reminiscent of the disorder induced phase control in such 
mixtures as discussed by one of us recently [T7| . In our 
case the phase control is not due to a disorder, but due 
to a control over the tunneling mechanism and the TRI 
breaking. 



II. BREAKING THE TIME-REVERSAL 
SYMMETRY 

A. One-dimensional optical lattice 

Let us begin, for simplicity, with a single particle in 
an one-dimensional (ID) optical lattice potential driven 
by a double harmonic perturbation. The Hamiltonian of 

the system reads 
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Ho = -^ + V{x)+Kixcos{ujt)+K2XCOs{2ujt + ip), (1) 

where V{x) = V{x + a) is an optical lattice potential 
with the lattice constant a, Ki^2 stand for strengths of 
the driving at a basic frequency w and its second har- 
monic, respectively. The Hamiltonian is time-periodic, 
i.e. Ho{t + 2n/uj) = Ho{t), and the Floquet theorem (isl - 
[20| guarantees that so called the Floquet Hamiltonian: 



v. = Hq — ihdt 



(2) 



is diagonalized by periodic functions. Eigenvalues of the 
H are referred to quasi-energies of the system analogous 



2 



to quasi-momenta in a solid state physics. They are de- 
fined modulo hu and it is sufficient to consider a single 
Floquet zone (an analog of the Brillouin zone). Peri- 
odic eigenfunctions of the Floquet Hamiltonian can be 
expanded in a Fourier series, i.e. in a basis where time 
t can be considered as an additional degree of freedom. 
Let us define basis vectors 



Ki . , , K2 . ,^ 
sin(wt) -f — — sm(2cji + ip) 



(\)j,m(x,t) = exp 

X exp (imLot) Wj (x) , 



(3) 



which fulfill 



27r/ti; 



{{(pj',m'\(l)j,m)) 



27T 



dt I dx ^rn'4'j,m 




(4) 



where m denotes a Fourier component and Wj{x) 



W(x 



is a Wannier function of the lowest energy 



band localized on the j-th lattice site. The first phase 
factor on the r.h.s. of Eq. ^ corresponds to the uni- 
tary transformation which allows us to switch from the 
length gauge to the velocity gauge using quantum optics 
language. The matrix of the Floquet Hamiltonian con- 
sists of diagonal m-blocks in the basis ([3]). The blocks 
are very weakly coupled among themselves provided the 
driving frequency cj is very high. Using the tight-binding 
approximation and taking into an account nearest neigh- 
bor tunneling only the diagonal blocks become 



(('/'j',m|^|'/'j,m)) — —Jc?i '^j'j'+l ^ '^off '^i'.j-l 



(5) 



where = {Wj\{p'^ /2m + V)\Wj) and the effective tun- 
neling amplitude 



JcS = J ^ J2k{si)Jk{s2)e 



ikip 



k——oo 



where 



aK, 

UJ 



(6) 



(7) 



are dimensionless strengths of the first (i = 1) and the 
second (i = 2) harmonic, the bare tunneling amplitude 
is J = —{Wj+i\{p^ /2m + V)\Wj) and Jn is the ordinary 
Bessel function, li hu ^ J the description of a single 
particle system may be restricted to a single diagonal 
block 



off 



(('/'i',o|'H|0j,o>> -^cff Sj'j 



+1 



(8) 



where the constant term Eq has been omitted. 

If there is only one harmonic present in ([T]) or the phase 
(ys = 0, the Floquet Hamiltonian is time-reversal invari- 
ant and H is represented by a real symmetric matrix in 



Figure 1: (Color online) The absolute value (top). A, and the 
complex phase ipj (bottom) of the effective tunneling ampli- 
tude, Joft/J — Aexp{i(pj), Eq. ([6]), for a double harmonic 
modulation of the optical lattice potential as a function of 
the dimensionless strength si of cu component for S2 ~ 1 [see 
(O], if = 0.2 (black solid lines) and S2 = 3, ip — 0.5 (red, 
dashed lines). 



a generic basis (2l|. Then, the effective tunneling ampli- 
tude ^ is real. Single harmonic driving has been used to 
change the interaction from ferromagnetic (positive Joff) 
to antiferromagnetic (negative Jeff) and to realize frus- 
trated magnetic phases (lll.[T2|. By breaking TRI we are 
able to realize nearly arbitrary complex values of the tun- 
neling amplitude Jeff = |Jeff|e*'^^ In Fig. [T] we present 
the absolute value | Jeff| and the phase ipj as a function 
of the parameter Sf. 

The eigenstates of the Hamiltonian ([5]) are the Bloch 
waves Ipj = e**''^'-' /a//Vs, where Ng is the number 
of lattice sites, with the dispersion relation E{k) = 
— 2| Jeff I cos(fca — ifj). Single harmonic driving allows for 
ifij = or TT and thus for the ground state with fc = or 
with k at the edge of the first Brillouin zone. The ground 
state of the system with broken TRI may correspond to 
any value of k. 

We have concentrated on a ID problem. However, a 
similar control of phases of the tunneling amplitudes can 
be also realized in higher dimensions. Indeed, the modu- 
lations applied to the lattice along orthogonal axes enable 
us to produce arbitrary tunnelings along the correspond- 
ing directions jTl| . In the following we will focus on a 2D 
triangular lattice. 



B. Two-dimensional triangular optical lattice 

The triangular optical lattice can be realized experi- 
mentally by means of three laser beams. Single harmonic 
modulations of the lattice along the two orthogonal di- 
rections allows one to control the sign of the tunneling 
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amplitudes of particles loaded in the lattice. This setup 
has been used in the Hamburg experiments that demon- 
strated the frustrated classical magnetism [l^]- With the 
help of the double harmonic modulations we are able to 
realize any phase of the tunneling amplitudes 



(9) 



see Fig.m Eigenstates of a single particle in such a lattice 
are Bloch waves with the dispersion relation 



^^(k) = —2\Ja\cos{kxa — Lpa) 

-2\Jf3\ |cos (ySky + k^'j I 
a 



^13 



(10) 



We induce a shift of the dispersion relation along the ky 
direction in the reciprocal space by changing the value of 
ipp (with the other parameters fixed). The modification 
of ipa alters the structure of the dispersion relation. It 
can reveal a doubly degenerate ground state for ipa = tt. 
The presence of such a degeneracy has been observed 
experimentally in a Bose system For example, for 

Ja = J 13 = ~|>^,3| the system in most experimental re- 
alizations chooses spontaneously one of the two ground 
states. With the double harmonic modulation breaking 
a TRI, the two degenerate minima for = tt can be 
moved arbitrarily along the ky direction with a change of 
ipp, see Fig. El 

In the Hamburg experiment a Bose-Einstein conden- 
sate (BEC) has been prepared in a triangular lattice. 
Although, in that case particle interactions are present, 
the ground state is still determined by the single par- 
ticle dispersion relation ([rU]) . Indeed, assuming an ho- 
mogeneous system (which is a good approximation of 
the experimental situation) the solution of the Gross- 
Pitaevskii equation has the chemical potential given by 
(jlb = E{\s.) -\- ubUb, where Ub characterizes the on- 
site particle interactions and ub is an average number 
of bosons per a lattice site. We would like to stress, that 
in the presence of the interactions, the restriction to a sin- 
gle block of the Floquet Hamiltonian like in ^ is valid 
provided huj ^ Ub [al- On the other hand, fuv must be 
much smaller than the energy separation between bands 
of the periodic lattice problem for the description limited 
to the lowest band to be valid. 



III. FERMIONS IN A TRIANGULAR LATTICE 

The frustrated classical magnetism in a triangular op- 
tical lattice has been demonstrated experimentally in a 
Bose system Within the mean field approximation 
Bose-Einstein condensate wavefunction is a Bloch wave 
with the wavevector corresponding to the minimum of the 
dispersion relation (|10p . For the antiferromagnetic inter- 
actions the system experiences the frustration, because 
the tendency of the wavefunction to change the phase by 




Figure 2: (Color online) Triangular Bravais lattice points 
(black circles) and amplitudes Ja^p corresponding to tunnel- 
ing from a lattice point to the nearest neighbors. 
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Figure 3: (Color online) Contour plots of the dispersion re- 
lation Eq. (|10[l for \Ja\ = \Jp\ and ipa = 1^13 — it/2 (a) and 
ipa — n and ipp — tt/4 (c); cool colors indicate regions around 
energy minima. In right panels directions of arrows indicate 
phases e''' ''' where k corresponds to a minimum of the dis- 
persion relation. Specifically ka = (f ' ^'^^ minimum 
in (a) and ka — ^+^, ^oy one of the two degenerate, 

non-equivalent minima in (c). Arrows in panel (b) and (d) 
are related to (a) and (c), respectively. 



TT. when we jump between neighboring sites, can not be 
reconciled with the triangular lattice geometry. 

Consider now a mixture of fermions in the different 
internal states (say spin up f and down J, states) with 
the attractive contact interactions in a 2D triangular 
optical lattice. We assume that the double harmonic 
modulations of the lattice allows us to adjust any phase 
of the complex tunneling amplitudes In the tight- 
binding approximation the Hamiltonian of the Fermi sys- 
tem reads 



Hp 



(»J> 
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Figure 4: (Color online) Modulus squared of the Fourier 
transform of the BCS pairing function |AkP, where Ak ~ 
Aie~^^''^^ /y/Ns, obtained numerically for a finite system 
of 60 X 60 lattice sites for |Ja| — |Js|, fa = n, Lp0 = n/4, 
U/\Ja\ = 2 and fi = 0. Panel (a) shows |Ak|^ corre- 
sponding to the ground state of the isolated Fermi system. 
Panel (b) presents similar results but for the Fermi system 
coupled to the Bose-Einstein condensate, with wavefunction 
ipi = ^/nse"''''''', where 'y^/n^/\Ja\ = 2.3. Note that the 
peak in (a) is located at ka = (0.00, 1.78) ~ ^0, ^1 while in 

(b) at ka = (2.09, 0.89) ^ qoa = (+'f, ^ 



-U 



(11) 



where the operator a^^ annihilates spin-up fermion at i- 
site, n.i^^ = dl^ai^ and similarly for spin-down fermions. 
The tunneling amplitude Jij = J*^ and it is equal or 
Ja, Eqs. ([5]), depending on a direction of the tunneling in 
the triangular lattice, see Fig. [51 The parameter U > 
characterizes the inter-species, on-site, attractive inter- 
actions and fi stands for the chemical potential of the 
Fermi system. 

The standard BCS approach (2^ leads to the effective 
Haniiltonian 



H 



F,cff 



(12) 



where the pairing function 

Aj = f7(a,;,-t-ai,^). 



(13) 



If the phases of the tunneling amplitudes ^ are zero the 
ground state of the system corresponds to constant pair- 
ing function A^ = const. However, the pairing function 
can acquire a non-trivial phase when the tunneling am- 
plitudes become complex. In order to find the ground 
state of the system let us look for the solutions of the 
Bogoliubov-de Gennes equations in the form 



(14) 



Wk(r») 










14 e"*'"' '"' 



where and Vk satisfy 



i;(k + ko) - 

A* 



A 




■f/k' 




■f/k" 


- ko) + /i 






= Ek 





(15) 

and |C/kP + |Vkp = 1. In Eqs. dH]), i;(k) is the dispersion 
relation pop while i?(k) = £'(k; ipa ~^ ~'^a, 
Solving (fl3)) we obtain the following eigenvalues 



-B(k + ko) - ^(k - ko) 



±fek 



(16) 



where 



lAP. 



(17) 

The excitation spectrum of the system, i.e. the upper 
branch ek.+; may become negative for some k. In such a 
case, the corresponding quasi-particles are present even 
at zero temperature. Therefore, at T = 0, the proper 
equation for A reads 



-E 



A 



2fek 



[1 - n- 



(18) 



where the Heaviside function 0(— ek,+ ) ensures that 
quasi-particles corresponding to negative energy spec- 
trum are also included {23| . Finally the desired pairing 
function becomes 



A, - C/^«k,+(r.X+(r.)[l-20(-ek,+) 



k 

„i2ko-ri^ 



(19) 



When we switch from (fa = '^p = to ifa = and 
ip/3 the minimum of the dispersion relation (jlOp is 

shifted from k = to k = ko = (^0, • In the ground 

state of the system fermions occupy energy levels starting 
from the new minimum up to the Fermi level. Thus, all 
fermions acquire quasi-momentum ko and consequently 
the pairing function gets the quasi-momentum 2ko. see 
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For |Jq,| = \Jp\ and (^q, = tt there are two non- 
equivalent, degenerate minima of Ei\<L). For example, 

for ipjj — 7r/4 they are located at k = (± 



■ 3a ' 20^3 

see Fig. In the ground state, fermions occupy en- 
ergy levels with quasi-momenta around both of the min- 
ima. The non-zero pairing function exists for differ- 
ent values of ko. However, for ko = ^0, 2a^z ) 
tain the lowest energy of the Hamiltonian p^ . A slight 
change of ko causes a rapid decrease of the energy gap 
in the excitation spectrum ([TS]). In Fig. we present 



Fourier transform |Ak| of the pairing function, where 
Ak = J2i Aie~*''''"'/-\/-^s, obtained numerically for a fi- 
nite system. For the chemical potential ^ ~ 0, that corre- 
sponds to half filling regime for non-interacting particles, 
\Jfi\/\Ja\ = 1, U/\Ja\ = 2 we obtain \A,\/\Ja\ = 0.111 at 
the center of the lattice which agrees with the analytical 
solution A/I Jq, I = 0.109 for infinite lattice. In Fig. we 
see that even for the lattice of 60 x 60 sites there exists 
a clearly resolved peak at k w 2ko = ^0, ^j^^ ■ 

It was discussed [l3|, Bose systems can simulate the 
frustrated classical magnetism. We show that the simi- 
lar phenomenon can be simulated in the Fermi system. 
Indeed, in triangular lattice with complex tunnelings, the 
phase of the complex pairing function is the one that is 
mapped onto orientation of classical spins. 



IV. BOSE-FERMI MIXTURE IN A 
TRIANGULAR LATTICE 

In this section we consider a situation when the 
fermions coexist with molecular dimers - pairs of spin- 
up and spin-down fermions. The dimers form a Bose- 
Einstein condensate. Such a mixture can be prepared 
by sweeping the system over a Feshbach resonance that 
creates a molecular BEC and leaves some fraction of un- 
bound, repulsively interacting fermions. Then crossing a 
second Feshbach resonance one is able to change the in- 
teractions between fermions from repulsive to attractive, 
turn unbound fermions into BCS pairs (l7| . The process 
does not affect the molecules at the same time. For this 
purpose Feshbach resonances at 202 G and 224 G for 
atoms [1^ seem to be quite suitable. We also assume the 
presence of a weak coupling that transforms dimers into 
unbound fermions and vice versa. It can be realized via 
the photo-dissociation and photo-association. For a large 
molecular BEC the weak coupling does not influence sig- 
nificantly the condensate wavefunction and therefore we 
neglect dynamics of the BEC. The system under our con- 
sideration can be reduced to the following Hamiltonian 



H — Hp + Hbf, 



(20) 



with 



Hbf = 7 X! {'^i "-ii^-it + i'i al-f 



^\ ' (21) 



where the BEC wavefunction ipi = ■y/nee*''''''"' is the 
ground state solution for bosons in the triangular lattice, 
i.e. qo corresponds to the minimum of the dispersion 
relation (|10p . For simplicity reasons but without loss 
of generality we choose the same dispersion relation for 
molecules and for fermions. In the system under con- 
sideration, the tunneling amplitudes for molecules in a 
shaken optical lattice depend on a molecular state pop- 
ulated in the photo-association process. The details of 
this process are not considered in the present paper. The 
coupling constant 7 characterizes transfer of dimers into 
unbound fermions and vice versa. We consider real 7 > 0. 

In the presence of the condensate of dimers the BCS 
effective Hamiltonian has to be supplemented with 
(EH), that is 



cff 



H 



F,cS 



H 



BF- 



(22) 



In the presence of bosons, if kg = qo/2, the simple 
analytic solution of the corresponding Bogoliubov- 
de Gennes equations exists. This solution need not cor- 
respond to the ground state of the system. However, 
we will see that for sufficiently strong coupling between 
bosons and fermions such a solution becomes the ground 
state solution. Employing with kp = qo/2 we obtain 
the following equation for A 



A 



^E^^^[l-2^(-^^+)]' (23) 



k 



where, in the present case, the excitation spectrum is 



ek,+ 



S(k + qo/2) -^(k- qo/2) 



+ fek, (24) 



with 



iE{k + qo/2) + ^(k - qo/2) - 2^^)= 



+ \A + -fy/W^f]^^^ , 

and the resulting pairing function 
A, = e*'i«""'A. 



(25) 



(26) 



Let us concentrate on the triangular lattice with 
I >^/3 1 / 1 -^ce I = 1; '/'a = TT and (pi3 = 7r/4 that corresponds 
to the dispersion relation plotted in Fig. The dis- 
persion relation reveals the two non-equivalent minima, 
but solution of the Gross-Pitaevskii equation for bosons 
chooses the Bloch wave with the quasi-momentum corre- 
sponding to one of the minima. The signatures of such 
a spontaneous symmetry breaking are observed exper- 
imentally We assume that Bose system chooses 

qo = , 20^3 } ^^'^ analyze its influence on the Fermi 

system. 
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We consider the system with /i = 0. If 7 = the 
Cooper pairs with the quasi-momentum qo do not ex- 
ist, i.e. A = is the only solution of (|23|) . If the 
coupling between bosons and fermions is present, but 
ly/riB /\Ja\ < 2.112, the system reveals gap-less superflu- 
idity [23|. The Cooper pairs with the quasi-momentum 
qo appear (A 7^ 0), but there is no energy gap in the ex- 
citation spectrum. The system possesses quasi-momenta 
k for which the excitation energies ek.+ < and con- 
sequently the corresponding quasi-particles are present 
even at T = 0. Concerning the ground state of the sys- 
tem, numerical solutions of the Bogoliubov-de Gennes 
equations are analyzed. It is found that increasing pa- 
rameter 7 causes a gradual enlargement of the peak at 
k = qo in the Fourier transform of the pairing function 

together with a reduction of the peak at k = ^0, 

(the solution in the absence of bosons considered in the 
previous section). For jy^nB /\Ja\ ~ 0.3 we observe a 

crossover, the peak at k = ^0, j becomes hardly vis- 
ible and the ground state starts to be well reproduced by 
the paring function (j26p . 

For jy^UB /\Ja\ > 2.112 the energy gap shows up, 
ek.+ > 0. There is no quasi-particle at zero-temperature 
and the pairing function (j26p is related to the ground 
state of the system. In Fig. HJd we show the Fourier 
transform of the pairing function obtained numerically 
for the triangular lattice of 60 x 60 site where the strong 
peak at k « qo is clearly visible. The pairing function at 
the center of the lattice is |Ai|/| Jq,| = 0.891 and the en- 
ergy gap in the excitation spectrum is 0.187| Jq|. Those 
numbers agree with the solutions for infinite system, i.e. 
A/I J„| = 0.891 and min(ek,+) 0.194| Ja|. 

Thus, we can describe the behaviour of the system in 
the following way. In a triangular optical lattice with 
complex tunnelings we are able to realize a BEC in the 
ground state with the wavevector located at an arbi- 
trary position in the reciprocal space. If the superfluid 
fermions are also present in the lattice and there is the 
sufficiently strong coupling between fermions and bosons. 



the phase of the BCS pairing function reflects the phase 
of the BEC wavefunction. 



V. CONCLUSIONS 

In summary, we have shown that the time-reversal 
symmetry breaking in an optical lattice potential allows 
us to realize complex tunneling amplitudes in the cor- 
responding tight-binding model. We have considered a 
simple scheme of the symmetry breaking by means of two 
harmonic modulations of the lattice, but the generaliza- 
tion to more complicated modulations is straightforward. 

We have studied a fermionic system as well as a Bose- 
Fermi mixture in a triangular lattice potential with com- 
plex tunnelings. In such a lattice the Bose system can 
simulate the frustrated classical magnetism [l2|. We have 
shown that this behaviour is similar for fermions where 
the pairing function acquires a complex phase. Assuming 
the presence of a coupling mechanism - an exchange of 
unbound fermions and bosonic molecules, we have shown 
that the complex phase of Bose wavefunction is mapped 
to fermions as reflected in the Fermi pairing function. 

While preparing the second, decompressed version 
of this manuscript we became aware of a very recent 
preprint where the authors consider similar idea for 
realization of complex tunneling amplitudes in Bose sys- 
tems both theoretically and experimentally. 
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